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SUMMAEY 

The variation of pressure distribution is calculated for a two- 
dimensional supersonic airfoil either experiencing a sudden angle— of— 
attack change or entering a sharp-edged gust. From, these pressure 
distributions the indicial lift functions applicable to unsteady 
lift problems are determined for the two cases. A close similarity 
is shown to exist between the calculated functions for varying free— 
stream Mach number and the corresponding functions in three- 
dimensional incompressible flow for varying aspect ratio. 

Be suits are presented which permit the determination of 
maximum increment in lift coefficient attained by an unrestrained 
airfoil during its flight through a gust. As an application of 
these results, the minimum altitude for safe flight through a 
specific gust is calculated for a particular supersonic wing of 
given strength and wing loading. 


INTRODUCTION 

The study of the unsteady lift of wings in an incompressible 
medium has been developed along two different lines. In refer- 
ence 1, E. T. Jones introduced the concept of indicial lift 
functions for wings of finite aspect ratio and, using as a basis 
the work of Wagner (reference 2) on the two-dimensional potential 
theory of airfoils in nonuniform motion, has shown how the 
calculation of lift under various conditions of motion can be 
effected. In reference 3, Theodorsen considered nonsteady motion 
in its relation to the general theory of aerodynamic instability 
and the determination of the aerodynamic forces on harmonically 
oscillating airfoils. This latter approach has been extended to 
include high-speed problems and in two recent papers Garrick and 
Eubinow (references k and 5) have given results on flutter and 
oscillating air— force calculations for wings in supersonic flow. 
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The present report employs the method of attaok Introduced hy 
Jones and considers the case of a two-dimensional airfoil moving 
supersonioally in an arbitrary manner, provided the assumptions of 
small perturbation theory are satisfied. The principal contribution 
lies in the determination of indicial pressure distributions which 
are readily calculated in supersonic motion and from which indicial 
lift, drag,, and pitching moments may be computed. From these 
results the indicial lift functions are calculated explicitly. 

The methods used to find the pressure distributions also afford 
considerable insight into the same problem for airfoils at subsonic 
speeds . 

As an application of the analysis, the results are applied to 
the special case of an unrestrained airfoil entering a sharp-edged 
gust. The resultant forces are found to be comparable in magni- 
tude, for Mach numbers in the neighborhood of 1.3# to those given 
in reference 1 for subsonic incompressible flow. 


SYMBOLS 

a speed of sound 

o chord length 

C} section lift coefficient 

C} a indicial lift coefficient for angle— of— attack change 
ct indicial lift coefficient for wing entering gust 

O 

g acceleration of gravity 

M free— stream Mach number 

m mass 

Ap difference in pressures between lower and upper surfaces of 

airfoil 

<1 (£p 0 v 0 2 ) 

s distance measured in half— chord lengths 
S area of wing 

p perturbation static pressure 
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t* t* t ime in seconds 

t transformed t ime variable (See equations (T)») 

T t/c 

u perturbation Telocity component in x direction 

V Q free-stream. velocity 

w perturbation velocity component in z direction 

w 0 z component of velocity of gust 

W weight of wing 

x 5 z Cartesian coordinates 

a angle of attack 

(i 2m/ pMSc 

p perturbation density 

p Q free— stream density 

a area over which surface integral, is evaluated 

<E perturbation velocity potential 

Subscripts 

0 free— stream conditions 

1 variable of integration 

u upper surface 

l lower surface 


AMLYSIS 

Derivation of Basic Differential Equation 

The pressure distribution over an airfoil in a compressible 
medium is obtainable from the solution of a boundary value problem 
associated with a particular second— order partial differential 
equation. The derivation of the linearized form of this equation, 
obtained under the assumptions of small perturbation theory, 
proceeds as follows : 
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Let u,w "be perturbation velocity components parallel, 
respectively, to the Cartesian axes x', z' and denote by 
p perturbation pressure, by p perturbation density, and 
by a the velocity of sound. Then if t' denotes time, 

V Q is the constant free -stream velocity, and p Q is the 
constant free-stream density, the linearized Eulerian equations 
are 


du du 

bt* + Y ° d?" 


dv 


+ Vr 


dv 

dx 1 


1 dp 

P o dx* 


1 dp 
P Q dz» 


CD 


The linearized equations of continuity and state are, respectively. 


dp 

dt 1 


- + Vo^- + P„(^_ ♦ %-) = 0 
• dx» °\dx» dz^ 


(2) 


p p 

a 2 
o 


(3) 


After the Introduction of the perturbation velocity potential $ 
integration of equations (l) yields 


dt + T d_£_ 
dt* 0 dx* 


— p + const. 
Po 


while equations (2) and (3) give the expression 


+ V 0 £_)— J2— 
>dt* dx’/ p 0 a 0 s 



W 


(5) 
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The combination of equations (4) and ( 5 ) leads to the desired 
partial differential equation 


(1_M2) d g| S + d 2 ® _ 2M df® 1_ 8 2 ® 

dx* 2 bz ,s ao 3x*St‘ a 0 2 8t* 2 


( 6 ) 


where M is the free— stream Mach number. 

Equation (6) can be reduced to the normalized form of the two- 
dimensional wave equation of mathematical physics by means of the 
transformation 


x = x* — lfe. 0 t* 

z = z * > 

t = at' 
o 


(T) 


In these variables the equation is written 

d 2 ® _ df® _ df® = 0 /g) 

3t 2 8x 2 8z 2 

In accordance with the assumptions underlying the derivation of 
equation (8), its application to problems in airfoil theory Is, 
of course, limited to cases where the induced velocities are small 
compared to the free— stream velocity and the effects of viscosity 
do not alter the results of the potential flow solution. 

The rectangular coordinate system associated with equation (6) 

Is fixed in the wing which Is, in turn. Immersed in a free stream 
of velocity Y 0 directed along the positive x-axis. The transforma- 
tion introduced in equation ( 7 ) fixes the x,z cordinate system 
in space so that the airfoil moves in the negative x direction and 
the free— stream velocity is zero. A distortion of the time axes 
is also involved so that the differential equation appears in 
canonical form. Equations (6) and (8) are, of course, well known 
in the theory of unsteady motion and in the study of sound waves . 

It is natural that these two fields of study should yield the same 
fundamental equations since, in small perturbation theory, the wing 
may be thought of as a distribution of acoustic radiators. 
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Solution For Given Boundary Conditions 

The boundary conditions which are to he satisfied have the Bams 
property as those encountered in steady— state thin-airfoil theory; 
that is, the prescribed data are given in the z = 0 plane. The 
particular problems with which this report^deals are those of 
finding pressure distributions over a flat plate. Thus, w will 
be specified over a portion of the z = 0 plane and, elsewhere 
throughout the plane, loading must be zero. 

Two boundary— value problems are to be considered: first, the 

case of an airfoil either starting from rest at a given angle of 
attack a or experiencing along the entire chord a change a in 
stream direction without a pitching motion, and second, the case 
of a constrained wing entering a sharp-edged gust with a vertical 
velocity w Q . In the former case the. motion is that of an airfoil 
suddenly sinking without rotation. These boundary conditions are 
more readily pictured with the aid of figure 1. Figure 1(a) shows 
the conditions which mustrbe satisfied in order to solve the angle— 
of— attack problem. The trace of the leading edge of the wing 
traverses the line x = — Mt, while the trailing edge lies on 
x = o — Mt where o is chord length. The region bounded by these 
lines and the line t = 0 is the region occupied by the airfoil 
as time passes. Since the exes are fixed and the airfoil moves in 
the negative x direction, the velocity at whioh the airfoil 
travels determines the inclination of the loci of the leading— and 
trailing— edge traces. Over the "area” occupied by the airfoil in 
the x,t plane, w must equal — Yoa and elsewhere no Jump in 
pressure can ocour. The gust problem (fig. 1(b)) does not differ 
essentially from the previous problem except that here the region 
over whioh the modification of w is effective is not entirely 
the region ocoupied by the airfoil but rather the region occupied 
simultaneously by the airfoil and the gust. Fixing, for convenience, 
the edge of the gust along the t— axis, this axis will form the right- 
hand boundary of the region over whioh w = — w D . 

Since the partial differential equation is linear and the 
solutions are therefore additive, these boundary conditions clearly 
should fit the following physical event: a wing of trapezoidal 

plan form. Indicated in figure 2, flies at a steady lift and angle 
of attack prior to t = 0; at t = 0 the wing either experiences 
a change in angle of attack a with no pitching motion or enters 
a sharp-edged gust of constant vertical velocity w Q , the gust 
extending from its edge to all negative values of x. The wing in 
each case is then restrained so that, relative to the original wind 
vector, the wing remains at an angle of attack a or, in the gust 
case, continues fixed at the same angle of attack. 
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The solution to similar ‘boundary— value problems has been 
discussed at length in reference 6. In that report the develop— 
nent was adapted to the case where equation (8) represented the 
steady— state equation for wing problems in three dimensions and 
the characteristic cones of the equation had the immediate physical, 
interpretation of Mach cones or infinitesimal shock disturbances. 
Green’s theorem was applied to solve the boundary— value problems 
involved and it was shown that the solution obtained could be 
interpreted as a surface distribution of sources and doublets. 
Finally* the difficulties arising in the discussion of the 
singularities on the characteristic cones and the integration of 
the supersonic doublet were overcome through the introduction of 
an integration technique which involved using the "finite part" 
of the given integrals. 

The methods and conclusions obtained in reference 6 can be 
adapted imme diately to the problems discussed herein for the 
mathematical reasoning remains almost identical. The physical 
interpretations of the two cases must* of course* be modified. 

Thus, the characteristic cones of the differential equations* 
traces of which are shown in figures 1 and 3, are no longer the well- 
known Mach cones; rather, they represent the distance to which a 
disturbance occurring at a point fixed by the apex of the cone 
will travel in the time t. Despite the fact that such physioal 
interpretations are undeniably useful in understanding and applying 
the results* the solution of the basic differential equation for 
the boundary values involved is quite independent of these material 
dissimilarities. It follows that the methods developed and discussed 
in detail in reference 6 can be applied directly to the given 
problems with only minor changes in notation. 

Deferring to figure 1, it is evident that for supersonic flight 
the air ahead of the wing is unaffected by the approach of the wing 
and* further* that the induced velocities on the upper surface of 
the airfoil are independent of the shape of the lower surface. 
Consequently the pressure distribution will be found on the upper 
surface* as if the airfoil section were symmetrical* and then* for 
the flat plate, the pressure distribution on the lower surface will 
be equal in magnitude and opposite in sign. As in reference 6* the 
solution to such a problem can be obtained from a distribution of 
sources. 

In the actual computation of the pressures over a section 
traveling at supersonio speeds* certain regions are conveniently 
defined. These regions depend on the relative slope of the traces 
of the leading and trailing edges and the trace of the character- 
istic cone in the x*t plane (fig. 3). The perturbation velocity 
potential is given by the form u la. 
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® 



dll dt^. 

(t— ti)2 — — z2 


(9) 


where a is the area in the x,t plane of the region occupied "by 
the wing section and hounded hy the trace of the forecone from the 
point x,z,t .. If p^ and p^ denote, respectively, pressures on 

the upper and lower surface of the wing 


AP 

<1 


p ? — P u 4/8® 8® \ 4 8® 

£ V 0 2 \8t' 0 8x’/ YqM 8t 


(10) 


Using equations (9) and (10) , direct calc ula tion showB that for a 
sudden angle— of— attack change, the following relations hold 

Region A (between lines x = — Mt, x = — t, and 1 = 0 — Mt) 


Ap _ 4a 
4 


(11a) 


It follows that in this region the steady— state Ackers t— type load 
distribution has been attained. 


Region B (between lines x * — t, x = t, and x 

Mx+t </ M 2 — 1 f % 


Ap 4a T 1 


arc cos + 

x+Mt 


xM 


■a 


+ arc sin 


= c - Mt) 

f)J < llb > 


Region C (between lines x = t, t = 0, 


and x = c — Mt) 


Ap _ 4a 
q. M 


(He) 


The result obtained for Region C is of particular interest since 
it holds for airfoils at subsonic as well as supersonic speeds. 
Moreover, the mechanics of the interaction between the airfoil and 
the fluid are such that other methods of derivation, furnishing 

added insight into the nature of the phenomenon, may be developed. 

Consider a flat plate of Infinite aspect ratio flying at a velocity 
Y 0 either greater than or less than the velocity of sound a Q in 
the undisturbed air. The airfoil is assumed to undergo a change in 
its motion at the time t = 0 so that subsequent to this t ime it 

has increased its angle of attack by the amount a. It follows that 

the sudden increment in lift can be calculated from a knowledge of 



MCA TIT Mo. 1621 


9 


the induced effects on the air produced "by an added vertical velocity 
of the plate equal to — Voa. As a result of this vertical motion two 
plane Rayleigh waves will emanate from the plate, a compression wave 
from the lower surface and an expansion wave from the upper surface. 
The velocities of the wave fronts are equal to a D while the 
induced velocities in the waves are equal in magnitude to Y 0 a. 

The lift on the wing can he determined in two ways : from 

impulse relations and from energy considerations. In the former 
case, assume that the forces per unit span on the upper and lower 
surfaces are f u and f j, respectively. After an elapse of time 
At the wave fronts have advanced a distance a^jAt and each 
includes a mass of air equal to pQcagAt per unit span. From 
Hewton’ s second law of motion 


( f i“ f u^ At = (p o 0a o At > (ar 0 *> 


and, converting to lift coefficient, 

= _ lta o a _ ^ 

° Z " ^P 0 V C " M 

Since the force is distributed uniformly along the chord this result 
is the equi valent of equation (lie) . 

In the development of the theory of plane waves of small 
amplitude (see, e.g., reference 7) Rayleigh and Lamb have shown 
that the energy in a wave is divided equally into kinetic and 
potential energy. Denoting kinetio energy by T, 

T = JJJ "* w 2 dxdydi 

where w is the perturbation velocity within the wave. Since the 
energy induced in the wave must result from work done on the plate, 
it can be seen that 
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or 

(f z -f u )At = 2p 0 Y 0 a ca^t 

This equation is in agreement with the one obtained, previously. 

As a consequence of equation (11c) it follows that the 
starting lift coefficient - of an airfoil is equal to Iki/m for both 
subsonic and supersonic flight. The magnitude of c i thus 
increases as M becomes smaller and for incompressible theory, 
where the velocity of sound is indefinitely large, must necessarily 
became infinite. This fact was known previously along with the 
understanding that the indicial lift function experiences an 
infinite discontinuity at t = 0. (See fig. 6.) For values of M 
other than zero the starting lift is finite and a continuous lift 
function results. 


For a sudden gust with vertical velocity w Q the following 
expressions can be found for the corresponding regions 


Region A 


fiefftpn.B 


Region C 


AP = w Q 
q v 0 «/m 2 -i 


Ap 

q 



Mx+t 

x+Mt 


AP 

— - 0 

q 


(12a) 


(12b) 


(12c) 


APPLICATIONS AND DISCUSSIONS 
Discussion of Load Distributions 

Figure ^ shows the variation of the loading on a section which, 
while traveling at supersonic speed, is suddenly deflected- to a new 
angle of attack. The loading varies according to equation (11) for 
each of the three regions A, B, and C of figure 3. At t = 0 the 
pressure is discontinuous. Jumping from, its original value. Just 
before the sudden deflection in angle of attack, to ifo/M Just 
after the deflection. Figure 3 shows, however, that the initial 
load distribution is modified over the forward portion of the section 
as time increases since regions A and B must be considered. 
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For t the wing lies entirely in region A so that 

Ap 4a 

<1 JWZL 

and the loading has attained a static value agreeing with the Ackeret 
type of distribution. The loading in region B varies between the two 
constant values of regions A and C, dipping below that of region C 
and having its minimum value at 2 = 0 . 

The loading produced upon entering a sharp-edged gust is pictured 
in figure 5 . The loading in region C is zero since that portion of 
the wing is unaware of the change in stream conditions. Over the 
forward portion of the airfoil (region A) the Ackeret type loading 
corresponding to the modified angle of attack is in evidence and for 

t >— 2— extends over the entire chord of the wing. In region B 
M— d. 

the loading experiences a reduction in magnitude from the value over 
the forward portion of the wing. 

The load distributions which have been developed were obtained 
for flight velocities in the supersonic regime. It is apparent, 
however, that the basic differential equation is not restricted to 
the case where M >1 and that the method of analysis affords a 
means whereby transient load effects can be studied for subsonic 
speeds. The essential difference between the latter problem and 
the results derived here lies in the relative position of leading— 
and trailing— edge traces in the x,t plane and the trace of the 
characteristic cone. Thus, for subsonic flight, the trace x = — t 
does not cut across the region occupied by the airfoil; whereas the 
cone stemming from the trailing-edge point c,0 does. A qualitative 
picture of the problem is obtained if the analogy between the non- 
steady two— dimensional case and three-dimensional wing theory Is 
used. The loading functions given in equation (11) are equivalent 
to loading existing on a swept— forward tip of a three-dimensional 
wing. Thus, in figure 1(a), x can represent distance measured 
spanwise, t can represent distance measured chordwise, and the 
shaded area can represent a portion of the plan form of the wing. 

Using this analogy, the loading which has been determined Is merely 
load distribution over the swept-forward tip of a wing with constant 
chord and supersonic leading edge. When the case of the airfoil 
section traveling at subsonic speeds Is to be considered, the problem 
becomes one of determining the loading over the swepi>-f orward tip of 
a wing with constant chord and subsonic leading edge . 
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Development of Indlclal Lift Functions 
Since section lift coefficient, cj,, is given by the expression 


C L 


l 

o J 0 q. 


dx 


the relations presented in equations (11) and (12) are sufficient 
for the determination of cj ■ , indicial lift coefficient for change 

in angle of attack, and CLg, indioial lift coefficient for an 

airfoil entering a gust. As a result of direct integration, the 
following results are obtained. 


First time interval 


tat, 

°La = TT 






cV, 


Second time interval 




1+M Mr-1 


( 13 a) 

(ItaO 




tat [* 1 /"jt 
* L M V2 


+ arc sin ~rr—) + - ■ ■■ ^ 


aro cos 


t+Mo-tM 2 


+ -j-*/ t 2 — (c— tM) 2 J 


(13b) 


W 0 t /it o-Mt\ kv Q Mc+t^—M^t 

c T » — { — + arc sin J + arc cos 

« \2 t / *7 0 «/Kl 0 

(Hft) 

tlmg.. laterra-Jr j^< t 
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Values of the lift functions are plotted in figure 6 as 
functions of s, the distance traveled "by the airfoil measured in 

2Mt 

half— chord lengths where s = — — . The curves shown were calculated 

for values of M equal to 1.2, 1.31* and- 1.46, since the asymptotic 
values of ct,^ and CL g for the three cases agree with the values 

given in reference 1 for the subsonic wing at aspeot ratios of 00 , 

6, and 3* respectively. No direct analogy, of course, can he made 
between the two cases. It is, however, worthy of note that the 
variations in the indicial functions for the superspnio case are of 
the same order of magnitude as those found in the finite— span 
incompressible case. 

From, a knowledge of the lift function resulting from a sudden 
uni t angle of attack, it is possible to express the lift corresponding 
to a given variable motion by considering the given motion as being 
composed of infinites ima l steps and summing the lifts corresponding 
to each step. Mathematically, the problem corresponds to the use 
of the indicial admittance in determining the current response for 
an electric network and leads to the so-called superposition -theorem 
which can be written in the form 

OL(V) - «(o) OI^(f) +£ OI^ (f-T.) JSLfVT- 

= — / (t *— t* ) a (t*) It* (15) 

The primes an the variables in this equation indicate that true time 
is used. 


Motion of Airfoil in Gust 


The results which have been obtained will now be applied to 
determine the forces on an unrestrained airfoil entering a gust. 
Since the motion of the wing is not prescribed it becomes necessary 
to equate the dynamical forces in order to relate the variables 
involved, neglecting pitching moment and using Newton’s second law 
of motion. 


m 


dw 

dt 7 


^ forces 


(16) 


where w is the vertical velocity of the wing, m is the mass of 
the wing, and the forces to be summed result from the lift on the 



14 


MCA TN Ho. 1621 


wing and the impressed force resulting from the action of the gust. 

By means of equation (15) , equation (16) can. be rewritten in the form 


+■ 1 

Avr A F* \r 

“dF* + 4S 5tr^ CI <X (t’-T*)a(T«)dr* = ^ CLg <f)qs 


Introducing a change of variables such that 


t * T* 

T = — a G , Tl = — a c 


(17) 


and setting 


— = a , 
V 

the equation becomes, finally. 


4 = 


2m 


PqMSc 


4 is + h S Q CI ^ (Mi) a (Ti)dTi = % (T) (l8) 

Since a(o) = o, ‘equation (18) can be integrated to give: 

m rji 

tt "vI CL e ( T i )dT i + ^ £ °Ia (T_Tl) a (Tl) dTl = 0 (19) 

which is an integral equation of the second hind with a variable 
upper limit. The solution to this equation can be obtained quite 
satisfactorily by means of Liouville’s method of successive 
substitution. * Using the relation c^ = and performing the 

proper manipulations gives: 

T 

«L ‘ °L g <*> “la, (^1) c Lg CWh 

* (Ti-Ts) “ Lg tel®2 — 


(20) 


Equation (20) is known to converge uniformly for T< 


\x*J M 2 — 1 


and in the applications of this report the maxi mum lift was always 
experienced in the region of convergence . 

The values of lift coefficient cjj determined from the solution 
of equation (20) are shown in figures 7(a), 7(b), and 7(c) for 
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various values of |i and for M = 1.2, 1.31, and 1.46. Figure 8 
shows the variation of the maximum lift coefficient attained plotted 
as a function of the density* parameter p. for the same three values 
of M. Figure 9 furnishes a comparison between the maximum lift- 
coefficient increment given in reference 1 for an aspect ratio of 6 
and the corresponding value calculated in the present report for 

M = 1.31. The results are plotted as functions of pj = -jSu to 

correspond with the density parameter used by Jones. The 
correspondence which was noted for the indicial lift functions with 
Mach number replacing wing aspect ratio is still in evidence. 


Forces Developed on Given Wing 

As an example of the uses to which the results Just obtained 
can be applied, consider an airplane with wing of plan form such 
as the one shown in figure 2. Assume a wing loading of 40 pounds 
per square foot, a chord length of 8 feet, and let it be specified 
that the wing is flying at a Mach number equal to 1.2 and that the 
wing is built to withstand forces producing accelerations between 
—3 and 5 times gravitational acceleration. It is proposed to find 
at what altitudes the wing may be subjected safely to a gust 
possessing a vertical velocity of 50 feet per second. 

If F denotes total force on the wing, 

F = W + L g (21) 

where W is wing weight and Lg is the total lift produced by the 
gust. If is the maximum increment in lift coefficient 

attained in a unit gust and A is the acceleration factor measured 
in multiples of g, then equation ( 21 ) is expressible in the form 


A = 1 + 


Since it is required that 


? § P ° T ° 


|A— l] = 4, it follows that 
< 6.4 


(acl) 


max 


PoVn 


( 22 ) 


Figure 10 shows the limit curve of (N]l) tot plotted as a function 
of flight altitude. From a knowledge of (i = 2W/p Q gMSc, however. 
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the value of aothally attained by the wing enter jng the 

gust can be calculated. Such values sire also included in figure 10 
and indicate that, under the given conditions, the wing should not 
fly at an altitude less than approximately 28,000 feet. 


Ames Aeronautical Laboratory, 

national Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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x=-Mt Position of leading 
edge after time t 


w =-V c ,<*■ in shaded 
section 


x--t \ 
trace of >. 
characteristi c 
cones 


(a) Supersonic wing after a sudden angle of attack 
change . 


x 



Figure /. - Boundary conditions for supersonic wings. 
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\^Trace of Mach cones 
\ from leading edge tip 

\ 

\ 
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Figure 2 Type of plan form studied in analysis. 




Figure 3.- Sketch indicating locations of regions A , B and 
C and time intervals used in analysis . 


mCA TIT No. 1621 


19 



(a) Time interval I , O < t/c < i/(M+i ) 


AP 

Q 

Regions 

4oc 

A 

B 

Im*-/ 

1 


v - 


leading edge trailing edge 

(b) Time interval II, l/(M J fl)<t/c<l/(M-l) 



(c) Time interval HI, i/(M-i)< t/c 

Figure 4.-Variations of Ap/q with chordwise station after sud- 
den change in angle of attack . Regions defined in figure J. 
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(a) Time interval I, O < t/c < I /(M+ 1) 



(b) Time interval n, f/(M+0 < t/c < f/(M-l) 


AP 

Q 

4w» 

v»4m*-i 



Region 

A 


leading edge 


trailing edge 


(c) Time interval Hit l/(M-l) < t/c 

Figure 5.-Variation of A p/q with chordwise station after enter- 
ing gust. Regions defined in figure 3. 
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Figure 7. - Variatiation of increment of lift coefficient during 
flight through unit sharp-edged gust. (a)M=l.ZO. 
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Figure 7. - Continued . (b) M= L3i 
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* Half -chords traveled \s 


Figure 7- Concluded, (c) M-L46 
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lift coefficient , -L-( Ac L ) max 


2.0 



Wing density parameter , jjl = 2 m/pScM 


Figure 8. - Maximum increment in lift coefficient attained during flight through unit 
sharp -edged gust for various values of free- stream Mach number, M. 
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Figure 9 . - Comparison of maximum increments in lift coef- 
ficients during flight through sharp-edged gust between 
subsonic wing and supersonic wing. 



Figure iO Safe altitudes for wing with loading of 40 !b/ft z 
flying at M=120 , based on 50 ft/sec. gust velocity 
and maximum load factors of 5 g and -5 g. 



